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Abstract. We discuss structure formation in cosmic string the- 
ories. Emphasis is laid on the properties of the peculiar velocity 
field, non-gaussian features and clusters of galaxies. It is found 
that the predicted peculiar velocity field is quiet (v < 150 km/s 
for L > 30/i _1 Mpc) for models with (n h < 0.2) which is 
in disagreement with the results from POTENT at large scales 
(> 30ft- -1 Mpc), consistent, however, with the analysis of the pe- 
culiar velocity field of spiral galaxies by Giovannelli et al. fT(| . 
Using the Press-Schechter formalism we calculate the abun- 
dances of X-ray clusters. It is found that in CDM models a 
mass fluctuation trg = 0.6. ..0.7 is needed with ttoh = 0.1. ..0.2 in 
order to explain the current data. 



1. Introduction 

An outstanding problem in cosmology is to explain the observed structures 
in the universe Two paradigms have been proposed: the inflationary 
cosmology and topological defects, which are well motivated by particle 
physics Jl3|, ^8|, (?], Both could, at least in principle, be tested with ob- 
servations of the microwave background anisotropies and with deep redshift 
surveys, because they make different predictions of the matter distribution 
and of the spectrum of the anisotropies of the cosmic microwave background 
radiation (CMBR). 

We discuss here the cosmic string theory, its prediction of the proper- 
ties of the peculiar velocity field, the non-gaussian features in cosmic string 
models and the implications of non-gaussianity on the abundances of clus- 
ters of galaxies. We do not discuss power spectra and CMBR anisotropies 
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Table 1. The four representative cosmological models, k is the cur- 
vature parameter, Q — 8itGp/(3Hq) is the matter density parameter, 
Ao = A/(3Hq) is the cosmological term, H is the Hubble parameter (in 
km s _1 Mpc _1 ). i?o and R eq are the scale factor at the present time and at 
matter-radiation equality, respectively. 



and refer instead to the contribution by Richard Battye (these proceedings 
and references therein). 

The paper is organized as follows: in section 2 we describe the cosmic 
string network evolution on large scales, which is relevant for structure 
formation. In section 3 we discuss the properties of the peculiar velocity 
field and in section 4 the non-gaussian features in cosmic string theory. 
Our conclusions are given in section 5. 

We take the opportunity to include some points which could not be 
covered in the talk. 

2. Network evolution 

A network of cosmic strings, which are line-like objects, will originate in 
a phase transition if the vacuum manifold is non-simply connected. To 
know how cosmic strings could form structures, we need to understand the 
evolution of such a network. Therefore we discuss this in more detail in 
this section, where we focus on the case of local cosmic strings. 

Our approach is based on the modified "one-scale" model by Martins 
& Shellard p5| . In this (phenomenological) model, the string network is 
characterized by a length scale L and the RMS velocity vrms of the strings. 
The length scale L is defined by 

Poo = W 

where is the energy density contained in the long strings (i.e. the strings 
with curvature radius larger than the Hubble radius) and fj, the energy per 
unit length on the string^ The model can not only be applied to cosmic 

1 fi is related to the energy scale of the phase transition where the string network 
was created. 
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Figure 1. Behaviour of L/Rh as a function of log(R/R(t eq )), where R(t eq ) 
is the scale factor at the time of matter-radiation equality and Rh is the 
(time-dependent) Hubble radius, adapted from pa|. 



strings but also to line defects in the condensed matter context. The reader 
is referred to the paper by Martins & Shellard, where further material and 
justifications of the equations are given fL5fl . We note that the equations 
are only exact in flat space-times. However, curvature effects can easily 
been incoporated. Our results presented below do not depend strongly on 
curvature effects. Curvature leads only to a faster release of energy of the 
network. Our results are in agreement with [Q. 

The equations, which determine the parameter L and wrmSj are given 
by (we set the speed of light c = 1): 



dt v 
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where k is a parameter related to the small scale structure on the strings 
and c is a parameter describing the energy loss due to loop formation. The 
expansion rate H(t) of the universe obeys the Friedmann equation: 




Figure 2. Behaviour of L/Rh as a function of log(Ro/R) for the models 
1 and 4. R is the scale factor at the present time and Rh is the Hubble 
radius, see [ p5[ . 



where p = p m attor+Pradiaton is sum of the total matter density and radiation 
density, A the cosmological constant, k the curvature parameter and R the 
scale factor. 

We solve the network equations for four cosmological models summa- 
rized in the table 1 J25| . The models include an open model, a flat model 
with a cosmological constant, a model with a loitering phase (see the arti- 
cle by van de Bruck and Priester, these proceedings) and the Einstcin-de 
Sitter model for comparison. The behaviour of the length scale L is plotted 
in Fig.l and Fig. 2. 

It is usually argued that the cosmic string network evolves according to 
the scaling solution, which states that the ratio of the length L to the Hub- 
ble radius Rn(t) = l/H(t) is constant: scaling solution L/Rh = const. We 
found that the scaling behaviour is always found in the radiation dominated 
epoch, but only reached in the matter dominated epoch in the Einstein-de 
Sitter model. In all other models, the scaling behaviour can not be found 
after the radiation dominated epoch.^ One can see, that even in the matter 
dominated epoch the scaling behaviour can not be found. This result can 

2 A epoch is called radiation dominated if the radiation density dominates all 
other terms in the Friedmann equation. The terms "matter dominated epoch" or "A— 
dominated" epoch are explained in the same manner. 
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easily be understood: As emphasized by Martins & Shellard, the transition 
from the scaling behaviour in the radiation dominated epoch to the scaling 
in the matter dominated epoch is long |lq |. In the low density models 1, 
3 and 4, the matter dominated epoch is too short for the string network 
to settle down towards the scaling behaviour, before the curvature term or 
the A-term dominates the expansion of the universe. 

Of course, in models with f^o ~ 0.9 the string network reaches scaling in 
the matter dominated epoch. However, such a model is indistinguishable 
from an Einstein-de Sitter model from the viewpoint of structure forma- 
tion. 

Very important for structure formation is the transition-regime between 
the radiation to the matter dominated epoch. We refer to the contribution 
by Richard Battye in these proceedings. 

In the cosmic string theory three mechanism for structure formation 
are possible: accretion on loops, wake formation due to fast moving strings 
and filamentary accretion onto slow moving strings, see e.g. 0, |8|. We 
us restrict here only on one of these possibilities: the formation of wakes, 
which is suggested by detailed simulations of the network evolution, see the 
discussion in |E8fl and references therein. 



3. The peculiar velocity field in cosmic string models 

Every model of structure formation is confronted with the observed proper- 
ties of the peculiar velocity field of the galaxies. An astonishing observed 
feature is that very large volumes (R > 20 Mpc) flow coherently through 
space (bulk flows) with large velocities (ftmik > 100 km/s, depending on 
the scale), suggesting the existence of structures on very large scales, see 
e.g. and references therein. 

Here, we calculate the bulk flows predicted from the cosmic string the- 
ory. In cosmic string theory, the peculiar velocity field of matter originates 
due to the conical space-time structure of a cosmic string. If a cosmic 
string passes a particle it will get a velocity kick in the direction swept out 
by the string The velocity kick is given by 

Wkick = -iirG^sVsf = 3.8(Gfj,) 6 (j s v s )f km/s, (5) 

with 



/ = l + m(l--)' ( 6 ) 



The term / corresponds to the small scale structure along the string, where 
ji is the effective mass per unit length along the string and T is the string 
tension. v s is the string velocity and j s is the corresponding Lorentz-factor. 
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On this basis it is easy to calculate the peculiar velocity field: the 
peculiar velocity v of a given volume is the sum of the influence of many 
strings, which could influence the volume: 

v = fkick / I -Aikj, (?) 

where Wkick is the value of the kick velocity due to a string, A% is the am- 
plification factor, which describes the amplification of the velocity kick due 
to gravitational instability and fcj is an unit vector pointing into the direc- 
tion of the area swept out by string i. The numbers of strings which could 
influence the volume are calculated with the modified one-scale model by 
Martins & Shellard |L5]]. The amplification factors can be calculated with 
the Zeldovich-approximation and we assume that the strings are uncorre- 
cted and therefore < ki -kj >= 8^ , where Sij is the Kronecker symbol j2f|. 



700 




v_rms/v_kick 



Figure 3. Probability function of the peculiar-velocity field u rm s (in units 
of the kick velocity i>kick) in the flat A-dominated model (model 4). Af- 
ter 50000 realisations the simulations produce an Maxwellian distribution, 
which shows the Gaussian character of each component of the peculiar 
velocity field, adapted from p(|. The dotted line is a gaussian curve to 
calculate the expected bulk flow (eq.10). 



We note that the peculiar velocity field can be described by a gaussian 
distribution (i.e. each component of v is a random variable with a gaussian 
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distribution). This is consistent with the results by other authors (e.g. 
[p^[). The probability distribution of the rms velocity v rms is plotted in 
Fig. 3 for the flat A-dominated model (model 4) and in Fig. 4 for the open 
model (model 3). For the other models we observe the same behaviour. 

2500 | , , , , , , 
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Figure 4. The same as in Fig. 3 but here for the open model. The solid line 
is a gaussian curve to calculate the expected bulk flow (eq.ll). 

We calculate v rms from the probability distribution by fitting a gaussian 
(see Fig. 3 and Fig. 4). In the models (Table 1) we obtain a peculiar velocity 
at a scale corresponding to L(t eq ) given by (the length scales are calculated 
with H = 100km/(s-Mpc)) §6|: 

Vpec{L e q ~ 70 Mpc) c l oscd (8) 

= (460±200)(G/x) 6 (7 s u s )/ km/s, 

V pec (L eq PS 1 Mpc) E ds (9) 
= (1740 ± 760)(G^) 6 ( ls v s )f km/s, 

Vpec(L e q ~ 10 Mpc) A: flat (10) 

= (280 ± 120)(G/i) 6 ( 7s u s )/ km/s, 
(L eq » 10 Mpc) 

open (ii) 

= (80±35)(G M ) 6 (7X)/Ws- 
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The length scale corresponding to the time t eq is set to be 0.17? (teg) 



L eq w 1.1— /iQ 2 Mpc 



(12) 
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Figure 5. Expected bulk flows Wbuik in cosmic string theory as a function 
of redshift in a flat, A-dominated model (model 4) and an open model 
(model 3). Here we assume (G^i)e(^ s v s )f « 1. The curves would shift 
if this quantity differs from unity (see text for details). One can see that 
cosmic strings in open or flat cosmological models predict a "quiet" peculiar 
velocity field at large scales. Note that in the closed model 1 this is not 
the case: the peculiar velocities are constant over a scale of < 10000 km/s 
and of the order 450 ± 200 km/s. 



Because the primordial spectrum of the density fluctuations in cosmic 
string theories is nearly of a Harrison-Zeldovich type, the peculiar velocity 
field is roughly proportional to L _1 (see Fig. 5). From this follows, for 
example in the open model (model 3) with fig — 0.1 at a scale of L = 50/i -1 
Mpc, a bulk velocity around 20 km/s, whereas in the flat A-dominated 
model (model 4) the velocity is around 70 km/s. This is inconsistent with 
the POTENT-results (see e.g. 0), as far as (Gfx)e('y s v s )f is of the order 
1, and with the Lauer-Postman result pH , which suggest bulk flows on 
large scales. These results are, however, consistent with the observed values 
obtained by Giovanelli et al. ||l6|| , which suggest a much more quiescent 
peculiar velocity field. The closed model fits the POTENT results well at 
each length scales. 
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The fact, that many cosmic strings influence the matter in the universe 
and the peculiar velocity field led us to ask the question if we can then esti- 
mate the density parameter from the comparison between the density field 
and the peculiar velocities. Such investigations are based on the equation 

V • v = -PH8, (13) 

with ft — f2[]' 6 , which holds in linear perturbation theory for an irrotational 
fluic^]. However, even if there is a relation between velocity and density 
field, eq.(^) is derived from the continuity equation and one can give 
examples in which an overestimated f2o could be derived from eq.dHf) (see 
the discussion for the case of the explosion scenario in ||). In the case of 
strings one can easily show that there is an connection between density and 
peculiar velocity field of the form 

V • v = -0 eS H5, (14) 

with 

PcS = ^ — t ■ (15) 

Here, 5ij is the density fluctuation created by the string i,j and /3 e ff,i,j 
is an "effective" value which connects the density and the velocity field of 
the string see [^6). In general, (3 c g should be different from Qq 6 . This 
let us to conclude that in general on can not obtain Qo from a compari- 
son of the density and the velocity field if cosmic strings where responsible 
for structure formation j2f|. However, important for this conclusion is the 
knowledge of the length scale at which non-gaussianity becomes impor- 
tant^]. One has to find this length scale with numerical simulations. 



4. Non gaussianity in cosmic string models 

Avelino et al. investigated this question for CDM and HDM models in 
detail They calculated the non-gaussian scale i? n g, below which the 
non-gaussian statistics of the density field is very important. For scales 
larger than R ng , the statistics of the density field can be described by a 
gaussian distribution to high accuracy. We discuss here the case of CDM. 

In order to find i? n g, Avelino et al. calcuated the skewness and kurtosis 
of the smoothed density field (i.e. of the probability density function) 0. 
Because the kurtosis and the skewness are zero for a gaussian distribution, 
these statistics give a hint for non-gaussianity if they are non-zero. Avelino 
et al. also calculated the genus curves of the density field, which is also a 
good indicator for the topology of the density field. 

3 Cosmic strings might generate vorticity also at large scales. However, these quickly 
decay by cosmic expansion (oc R~ 2 ) so that we can neglect it at the current epoch. 

4 I am grateful to Pedro Avelino for pointing out this to me. 
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Figure 6. The tails of the effective probability function fP(y) for gaussian 
and cosmic string models on the scale i? ng . The cosmic string data are 
taken from [Q. y = 5/aR, where S is the density contrast and or is the 
mass fluctuation at the scale R n g, see the text for details. 



As a result they found R ng = 1.5(r2o/i 2 ) - 1 Mpc. For models with Slo ~ 1 
this length scale is very small compared to 50 Mpc at which large scale flows 
are observed. Therefore, for these models the density and peculiar velocity 
field are related by equation (jp]). In fact, only for models with floh 2 < 
0.1 the scale is large enough that a significant departure from eq.([l3|) is 
important (for example model 1). Then eq.(|l4|) holds. Cosmic strings in 
such very low density models might be ruled out by the correlation analysis 
in Fourier space, see Stirling & Peacock pl[ . They found from the combined 
QDOT and 1.2-Jy IRAS galaxy survey that the modes with wavelength 
A > 30/i _1 Mpc are in good agreement with a gaussian distribution. While 
this is a provisional conclusion, with future observations this test will place 
limits on the degree of non-gaussianity in the density fluctuation field on 
small and large scales. For another discussion on non-gaussianity see 

In low-density models with floh — 0.1. ..0.2 the length scale R ng is about 
10 Mpc. From this length scale clusters collapse. Therefore, different pre- 
dictions of the number density of cluster of galaxies are expected. In fact, 
clusters are used to test cosmological models, see e.g. || and references 
therein. Here we use the Press-Schechter method to caluclate the abun- 
dances of clusters in cosmic string theory. Because of non-gaussianity we 
have to use a modified approach, see e.g. |J. 
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We consider the smoothed density field 8. The probability density func- 
tion p(S) can be described by a rescaled function P(y) withy = S/aR, where 
<tr is the mass fluctuation at the scale R. We assume p(5) = P(y)/aR which 
should hold in the theories considered here and on the relevant length scales. 
The probability density function P(y) is shown for cosmic string models 
and gaussian models in Fig. 6. 

In the Press-Schechter formalism the number density n(R) of collapsed 
objects (collapsed from a radius R) is given by 



n{R)dR 



AttR 3 



fP(y)dy 



dR, (16) 



where / = 1/ J °° P(y)dy and y c = 5 c /cjr. Here 5 C is the critical density 
contrast for collapse. We assume the spherical collapse model with 5 C — 
1.68, because this quantity depends only very weakly on the cosmological 
model [ Ell . From the power spectrum we calculated <tr. It can be shown 
that (Jr oc R~ a , where a ~ 1. Our results for the range 1 keV to 10 
keV depend not very strongly on the slope of the power spectrum if a = 
0.75. ..1.25, which was also observed by others 111). We give details of our 



calculations elsewhere |27 



Because the length scale R from which a cluster collapsed is not an 
observable, we have to relate R with the cluster mass M or the X-ray 
temperature of the cluster gas. Here we calculate the temperature function 
and assume a spherical cluster collapse. The resulting cluster is described 
by an isothermal sphere with a radius of virialisation equal to half its 
maximum expansion radius (e.g. ]3C|]). Then, for a hydrogen mass fraction 
of 0.76, 



kT = 8.6 keV fi^ 3 ( —4 ^ 2 ( JO V ' 



ro/j-iMpcy \n(zf)J 



Ac 
178 



1/3 

(17) 



where z/ is the redshift of collapse, which we calculate with the spherical 
collapse model. We calculate the number density n(> kT) of clusters which 
X-ray temperatures greater than kT (temperature function). 

Because high density peaks are much more common in cosmic string 
models as in gaussian models, for a given density parameter, cosmological 
constant, Hubble parameter and mass fluctuation at 8 h~ l Mpc, cosmic 
strings predict more clusters with high X-ray temperatures as compared 
to gaussian models. This can be seen in Fig. 7. If we change as the tem- 
perature function changes (Fig. 8). Also shown are the data by Henry & 
Arnaud (1991) (the vertical bars with the range of uncertainty). We found 
that ag should lie in the range 0.6 — 0.7 to fit the data (for a flat model 
with fio = 0.3). The reader is referred to our forthcoming work, where more 
data will be used and where the errors, which could enter the calculations, 
are discussed in detail p7| . 
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Figure 7. Number densities of clusters of galaxies as a function of the 
temperature kT in a gaussian and a cosmic string model. The cosmological 
parameters are Qq = 0.3, A = 0.7, h — 0.7 and as = 0.75. 



Similar to gaussian models we find that the temperature function (and 
also the mass function) are good discriminators for non-gaussian models. 
We refer to our future work, where these functions and cluster evolution 
will be discussed in some depth f27|j . 

We are not able to discuss cluster abundances in the baryonic model 
1 of Tab.l because we don't know the probability function on the scale 
relevant for cluster formation. 

In conclusion we point out, that several uncertainties enter the analyis: 

• We concentrate here only on the wakes. How important are cosmic 
string loops for structure formation? Is filament formation due to 
slow moving strings also important? If these mechanism turn out to 
be important, the probability density function of the density field will 
certainly change. 

• How does the cluster gas evolve and what is its impact on the temper- 
ature function (e.g. |22))? Here, an detailed analysis of the formation 
and evolution of wakes seems necessary, see [^3[ [29| . 
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Figure 8. Number density of clusters of galaxies in cosmic string models 
with ilo = 0.3, Ao = 0.7, h = 0.7 and different as as a function of the X-ray 
temperature. Shown are also the temperatur function according to Henry 
& Arnaud @. 



• How good is the spherical collapse model? What is the typical for- 
mation redshift of a cluster? 

• Is the Press-Schechter formalism good enough do describe the statis- 
tics of high-density peaks in non-gaussian models? 

At the present time we are far from conclusive answers for the questions 
above. 

The case for HDM models + cosmic strings is not discussed here. But 
the simulations suggest that on cluster scales the probability density func- 
tion is gaussian in these models and cosmic strings may have the same 
problems as in "usual" gaussian models Qf]. Note, however, that it is im- 
portant to know which kind of perturbations are important in cosmic string 
theory: wakes, filaments and/or loops are possible, for a discussion see e.g. 

S0- 

5 However, see [ po| , for a discussion of the formation of high redshift objects in cosmic 
strings + HDM. 
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5. Discussion 



We have discussed some aspects of the cosmic string theory of structure for- 
mation. Our results suggest that low density models are very encouraging 
for cosmic string theory. Problems might arise if the peculiar velocity field 
is observed to be of the order of 300 km/s at large scales (cz > 3000 km/s). 
Then only models with fl^h < 0.1 and a cosmological constant could pro- 
duce such large flows at these scales (if {G^)^ s v s )f w 1). Such models 
predict non-gaussian features at these scales and could be constrained with 
topological statistics, see e.g. jl8) or with correlation analysis in Fourier 
space, see e.g. [ gjj . 

Further observational efforts will help to decide if cosmic strings can be 
reponsible for structure formation. In particular if the peculiar velocity field 
is better constrained from observations at large scales and if we understand 
cluster formation and evolution more thoroughly, we should able to test 
these theories apart from CMB observations and matter distribution. 
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